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Abstract 

The spin-lattice relaxation time tsl is a key quantity both for the dynam- 
ical response of ferromagnets excited by laser pulses and as the speed limit 
of magneto-optical recording. Extending the theory for the electron param- 
agnetic resonance of magnetic impurities to spin-lattice relaxation in ferro- 
magnetic rare earths we calculate tsl for Gd and find a value of 48 ps in 
very good agreement with time-resolved spin-polarized photoemission exper- 
iments. We argue that the time scale for t$l in metals is essentially given by 

the spin-orbit induced magnetocrystalline anisotropy energy. 
78.47.+p,75.50.Cc,79.20.Ds,75.30.Gw 
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I. INTRODUCTION 

The spin-lattice-relaxation time t$l is a sensitive fingerprint for the strength of the dynam- 
ical coupling between the spin system and the lattice. This time is therefore of interest for 
the long-time spin response of magnetic materials upon pulse laser lattice excitation. In a 
ferromagnetic solid, this time is required to establish a new equilibrium magnetization after 
a sudden change of the lattice temperature. Thus, t$l is a key quantity for magneto-optical 
recording, since it determines the maximum speed for magneto-optical Curie-point writingS. 
In their pioneering work on Gd, Vaterlaus et a/.i were the first to measure tsl in real-time 
using time-resolved spin-polarized photoemission. This experiment was performed with the 
pump and probe technique applying strong 10 ns laser heating pulses followed by 60 ps weak 
probe pulses with variable delay and yielded the result 

r SL = (100 ± 80) ps = 5 - 50GHz . 

This corresponds to a gain by two or more orders of magnitude in speed compared to the 
present state of the art of data processing, which is still rapidly improving. 

Up to now there exists no calculation or theoretical explanation of this result. Thus, it is 
the goal of this paper to provide a theoretical approach to spin-lattice relaxation in metallic 
rare-earth ferromagnets, which is based on rate equations and an electronic model structure. 
We present a theory which, despite of its simplicity, exhibits already the microscopic features 
of spin-lattice relaxation in these materials. 

The spin-lattice relaxation time t$l describes the time required by the spins to reach 
thermal equilibrium with the lattice. The lattice then operates as a heat bath if one neglects 
the "phonon bottleneck" thus assuming perfect coupling to the external environment via the 
phonons. Hereby the originally cold spins are flipped by the phonons, and spin and phonon 
systems approach a common thermal equilibrium. This is microscopically accomplished as 
follows: The spins couple to the anisotropic fluctuations of the crystal fields produced by the 
phonons. This coupling is mediated by spin-orbit interaction. During this process neither a 
modification of the geometrical structure nor a change of the magnetic phase (long range 
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order) has to take place. 

A typical scenario of the processes leading to spin-lattice relaxation is a four-step process: 
(i) The laser beam hits the sample and creates electron-hole pair excitations within 10~ 15 
sec. (ii) The electronic system equilibrates at elevated temperatures by electron-electron 
interactions within 10 fs. Note, the lattice is not yet involved. The spin and charge dynamics 
at the elevated temperature may already lead to the breakdown of magnetic long range order 
in the case of intense fs laser pulses but not for ns heating, since the electronic system is 
always close to equilibrium for ns photoemission. (iii) The equilibrated electronic excitations 
decay via phonon cascades within 10~ 13 . . . 10~ 12 s and heat up the phonon system, i. 
e. the lattice, (iv) The phonons and the spin system reach their common equilibrium 
within the spin lattice relaxation time tsl of 10~ 10 s. This is also the time which allows 
for the recovery of magnetism in many cases, since the electronic equilibrium temperature 
after step (ii) might be much larger than the Curie temperature, whereas the common 
equilibrium temperature reached after step (iv) for spins and phonon system is usually 
much lower than the electronic equilibrium temperature and may also be smaller than the 
Curie temperature. The characteristic interactions of these four processes taking place on 
distinct time scales are: (i) p • A, where p is the crystal momentum of the electrons and 
A is the vector potential of the laser photons, (ii) electron-electron Coulomb interaction 
leading to dynamical charge and spin fluctuations, (iii) electron-phonon interaction, (iv) 
phonon-magnon interaction caused by spin-orbit interaction which we will approximate by 
the static magnetocrystalline anisotropy energy (see below). 

The experiment by Vaterlaus et al. was the first to measure the time evolution of the 
magnetic nonequilibrium state on a picosecond time scale. Therefore, in this paper, we 
exclusively address the long-time (ps to ns) response via the lattice to compare with the 
above experiment. 

To support the above scenario of a laser pulse causing on the ps timescale and at not 
too low temperatures mainly the heating up of phonons, we compare the specific heat of 
phonons, spins, and electrons! and find the following: The spins start to dominate the 
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phonons at temperatures 



T < T « 0.1 6 D (1) 

for fields of about 1 Tesla (in the case of paramagnetic impurities). The electrons start to 
dominate the phonons at temperatures 

T < T « 0.01 B D . (2) 

Typical Debye temperatures Qd in ferromagnets are 420 K (Fe), 385 K (Co), 375 K (Ni), 
152 K (Gd), 186 K (Dy). (For details of these estimates see appendix A). This crude esti- 
mate shows already, that phonons are dominant at sufficiently long time scales and not too 
low temperatures (the experiment has been performed typically at temperatures between 
30 and 300 K), whereas spins (corresponding to Coulomb-correlated electrons) and finally 
electrons (single-particle excitations) are going to take over for lower temperatures but also 
for shorter times (100 fs and shorter). Thus, it appears reasonable to focus on the phonons, 
since the experiment has been done on the ps to ns time scale and at low temperatures. 
However, it becomes immediately obvious that very interesting dynamical properties of the 
electrons and spins are to be expected in faster (fs) pump and probe experiments which 
will definitely be available in the near future. However, for this time window, the notion 
of spin- lattice relaxation makes no sense any more, since electrons rather than phonons are 
involved. 

II. MICROSCOPIC CALCULATION OF THE SPIN-LATTICE RELAXATION 
TIME T SL 

To calculate the spin-lattice relaxation time t$l we start from the theoretical approaches suc- 
cessfully applied to electron spin resonance (ESR) more than three decades ago for magnetic 
impurities embedded in a nonmagnetic host lattice and adapt this treatment to the solid 
combining phenomenological nonequilibrium thermodynamics (kinetic theory) and micro- 
scopic equilibrium theory. Three processes (all involving phonons) contribute to spin-lattice 
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relaxation: The direct process (Fig. 1 (a), see appendix B), the Orbach processo (Fig. 1 
(b), see appendix B), both of them being relevant only at very low temperatures, and the 
Raman process (Fig. 1 (c)) which we consider here: This process consists of a spin-flip, 
the absorption of a phonon of frequency u, and of the emission of a phonon of frequency 
uj + The longitudinal relaxation rate T\ in this case is independent of the magnetic 
field! and is given by 

— ~ T 7 . . . T 9 . (3) 

The Raman process is a two-phonon process of higher order which essentially uses the 
complete phonon spectrum. This process dominates the Orbach process (and thus also the 
direct process, see appendix B) for 

|>e D , (4) 

where Ai is the crystal field splitting and ks is Boltzmann's constant. Nickel, for example, 
has 

^ « 688K > <d D » 375K . (5) 

Thus, for not too low temperatures, the Raman-process is dominant for the spin-lattice 
relaxation rate. 

Therefore, in view of the experimental conditions, it appears justified to focus on Raman 
determined spin-lattice relaxation in the solid which should be valid at intermediate lattice 
temperatures and ps time scales. The temperature range of validity forms probably the 
best compromise between too large temperatures where the lattice becomes unstable (above 
the melting point) or magnetism breaks down (above the Curie temperature) and too low 
temperatures where direct and Orbach processes determine the phonon induced relaxation 
or the phonons become frozen. Besides, the Raman process is independent of the magnetic 
field. 

Note that purely electronic mechanisms such as spin fluctuations in strongly correlated 
electronic systems mediated by nuclear spin-flips (for energy and angular momentum con- 



servation) via hyperfine interaction require even longer time scales and are unimportant in 
this context since they do not involve the lattice. 

To calculate now Raman-induced spin-lattice relaxation in ferromagnetic rare earth solids 
we start from the theory for spin-lattice relaxation in magnetic impurities^. First we consider 
the number of phonons in the volume V and energy interval [5,5 + d5] 

/x . 3V6''d6 

where v s is the speed of sound in the material (e. g. Gd). The thermal occupation is given 
by the Bose factor: 

Po(S) = ~^g~ ■ (7) 
e k B T — 1 

For the interaction, the usual crystal field expansion up to second order in terms of the 
randomly fluctuating strains is used 

K « ei£2$>n> (8) 

ran 

since the Raman effect is of second order (see Fig. 1(c)). The transition probability from 
state 1 6 > to \a > is then given by 

w b ^ a = J — \< b,po(S 1 ),p (S 2 ) | Hi | a,p (5i) - l,p (5 2 ) + 1 >| 2 p(S 2 )p(5 1 )dS 1 . (9) 

Including the processes of stimulated emission, absorption, and spontaneous emission the 
rate equation for the change of the occupation numbers of the levels \b > and \a > is given 
by (p is the mass density of the solid) 

N b = -N b w b ^ a + N a w a _ b = -N a = K[-N b po(5) - N b + N a p (5)}. (10) 

Using eqs. (8) and (9) this leads to 

9 V \< n I v m I h >l 2 r 
^ = J 2 3 |7 n J ~ / [NMh)\po(5i) + 1] - A^o(5i)[p Q (<y + 1]]^. (11) 

Here, it has been used that the square of the matrix elements of the strains e assumes the 
value 



SjMS) + 1] 
2Mvl 



(12) 



where M is the crystal mass. Using the plausible assumptions 

5 < k B T, 5 < «5i (13) 

and the abbreviations 

n = N a — N b , N = N a + N b and n = iVtanh( — — ) (14) 
yields the kinetic equation of spin-lattice relaxation 

h = (n — no) ■ (15) 

T~SL, Raman 

The microscopic calculation of the spin-lattice relaxation rate (which is the kinetic coefficient 
of the rate equation) gives then the result 



1 9Em« \<a\v™\b >| 2 [ k B o Sfe^ddi 



TSL,Raman 8 p 1 ''TT 3 h' 'v] 



(16) 



Using our previous estimate for the magnetocrystalline anisotropy, which is discussed in 
some detail in appendix C, 

]T |< a | < | b >| 2 = | E anisotropy | 2 = | 735/ieU | 2 , (17) 

mn 

this microscopic theory finally yields for the spin-lattice relaxation time in Gd a value of 

T S L,Raman = 48pS . (18) 

This result is in excellent agreement with the experimental value of (100 ± 80) ps. The main 
issue here is that obviously the energy scale for spin lattice relaxation is set by the magne- 
tocrystalline anisotropy energy, which is of the order of 100 fieV - 1 meV at surfaces, in thin 
magnetic films or in hexagonal bulk crystals, rather than by the Curie temperature or by 
spin-orbit coupling or by electron-phonon interaction (all being of the order of 30 - 50 meV). 
This energy scale comes into play, since spin-lattice relaxation orginates from the coupling 



of the spins to the anisotropic crystal field fluctuations resulting from the phonons. These 
fluctuations flip the spins to accomodate their thermal occupation to the lattice tempera- 
ture (or to a common equilibrium spin- lattice temperature). Although magnetocrystalline 
anisotropy results from spin-orbit coupling, its energy scale is typically smaller at interfaces 
or in the bulk of noncubic threedimensional solids by a factor of 100, since spin-orbit cou- 
pling enters to second order (see appendix C). In cubic bulk crystals the leading terms are of 
fourth order thus resulting in a reduction factor of 10000. This argument holds for both (i) 
the level shifts induced by spin-orbit coupling and (ii) the occurrence and lifting of degenera- 
cies at the Fermi energy within a small portion of the Brillouin zone!. Our argumentation is 
still valid even for the particular case of Gd, where the localized /-shell carries most of the 
magnetic moment while the conduction electrons are responsible for the metallicity, since 
the anisotropy of the magnetic moments involves the coupling of localized and conduction 
electrons. The same holds for the spin-lattice relaxation. 

III. CONCLUSIONS 

In this work we presented a microscopic theory for the spin- lattice relaxation time tsl m 
the metallic rare earth ferromagnet Gd and found a value of 48 ps in remarkably good 
agreement with experiment. Although our theoretical estimate neglects all detailed features 
of electronic structure, phonon density of states, electronic correlations, effects of electronic 
temperature, and the detailed form of the transition matrix elements it already yields the 
correct value of t$l- Moreover, our theory clearly demonstrates the important relationship 
between the static magnetocrystalline anisotropy energy and the dynamic quantity tsli 
which is essential for magneto-optic recording velocities. Furthermore, our theory yields a 
good starting point for a detailed electronic and nonequilibrium response theory of spin- 
lattice relaxation in rare earth and transition metals (involving phonon-magnon coupling, 
see appendix D). Thus it could overcome the restriction of previous ESR theories to localized 
magnetic impurity spins (e. g. in insulating garnets). For thin films and multilayers, it is of 
particular importance to calculate the thickness dependence of this dynamical quantity, thus 
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checking the range of validity of the relation of tsl to magnetocrystalline anisotropy and 
to the linear and nonlinear magneto-optical Kerr-effects, which involve the complementary 
non-spin-flip effects of spin-orbit interaction. It is a considerable theoretical challenge to 
investigate phonon-magnon coupling in realistic itinerant systems. From the experimental 
side, additional thickness dependent and spectroscopic pump and probe laser experiments 
as well as measurements of the ferromagnetic resonance (FMR yielding collective spin-flip 
frequencies) are required to tackle the important and complex problem of spin-lattice re- 
laxation in metallic ferromagnetic thin film media and to bridge the gap between magnetic 
resonance experiments in the frequency domain and optical real-time measurements. In par- 
ticular, it will be interesting to study the temperature dependence of t$l, thus discussing 
also low-temperature contributions to the relaxation originating from direct or Orbach pro- 
ceses. 

Besides, it is of considerable interest to search for faster spin-switching mechanisms using 
intense fs laser pulses which may directly lead to a breakdown of magnetism via electron- 
electron correlations and may therefore bypass the lattice thus reducing lattice heating. It 
is to be expected that more interesting results will be found on the femtosecond time scale 
which is now also accessible using Ti-sapphire lasers. Upon intense laser excitation, the mag- 
netic state may break down already within some fs without the influence of the lattice and 
it is recovered within tsl which involves coupling of the spins to the lattice via anisotropic 
crystal field fluctuations. In this case, the spins are cooled by the lattice rather than heated 
as in the experiment by Vaterlauset ai, which requires a theoretical explanation. These 
time scales should be optically accessible in metallic thin film media in the near future. 

ACKNOWLEDGMENT 

We gratefully acknowledge stimulating discussions with Prof. K. Baberschke, Dr. P. J. 
Jensen, and Dr. F. Meier. 



9 



APPENDIX A: SPECIFIC HEAT OF PHONONS, SPINS, AND ELECTRONS 

In this appendix, we compare the specific heat of phonons, spins, and electrons! in order to 
support our approximation that, at not too low temperatures, a 10 ns laser pulse heats up 
mainly phonons. 

The low-temperature specific heat of phonons at constant volume is given within the 
Debye model by 

127r^ T T ^ 

c v = ~^-rik B (—) 3 = 234— nk B . (19) 

O <&£) C>i3 

Here, Qu denotes the Debye temperature, k B is Boltzmann's constant, and n is the number 

of lattice sites per unit volume. 

The specific heat of the spins is given by 

" = 57 fe ' J < J + 1 ><W ) "' (20) 

where J is the total angular momentum and g is the gyromagnetic ratio. 
The specific heat of the electrons is 

( dn ir 2 2 n 2 k B T 
cv = (gy) = —k B Tp{e F ) = y— nk B . (21) 

The spins start to dominate the phonons at temperatures T 

r<T. = (^)*(^)«e D »aiex, (22) 

for fields of about 1 Tesla and iVj paramagnetic ions. The electrons start to dominate the 
phonons at Temperatures T 

T < T = 0.145(4^)^e D « 0.01 Q D , (23) 

J-F 

where Z is the nominal valence. Typical Debye temperatures in ferromagnets are ®d = 
420 K (Fe), 385 K (Co), 375 K (Ni), 152 K (Gd), 186 K (Dy). This crude estimate shows 
already, that, for ns laser pulses, phonons are dominant at not too low temperatures. 
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APPENDIX B: DIRECT AND ORBACH PROCESSES 

In this appendix, we discuss the direct and Orbach relaxation processes which may dominate 
the Raman contributions to spin-lattice relaxation only at very low temperatures, 
a) Direct process (Fig. 1): 

This process consists of a spin-flip and the emission of a phonon of frequency u . The 
longitudinal relaxation rate is proportional to the temperature 

The rate is also proportional to the number of phonons within a narrow interval S at the 
extreme low-frequency end of the phonon spectrum 

MS) « (25) 

Thus, the direct process is of importance only for very low temperatures (T <C ©d)? where 
the other processes become negligible. The relaxation rate for the direct process depends on 
the magnetic field (for magnetic impurities). It is proportional to H 4 for Kramers-doublets 
and proportional to if 2 for non-Kramers-doublets. 



b) Orbach process (Fig. 2): 

This processS consists of a spin-flip, the absorption of a phonon of frequency 4, and the 
emission of a phonon of frequency j- + u . The longitudinal relaxation rate for the Orbach 
process is approximately given by 

1 A i 

— ~ e k B T . (26) 
T\ 

This rate corresponds to two high frequency cascades and is proportional to the number of 
phonons in a narrow band at 

Ai w crystal — field splitting 

and does not depend on the magnetic field. The Orbach process becomes important if the 
relation holds: 
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(27) 



At higher temperatures, however, such as in the experiment by Vaterlaus et ai, the Raman 
process should dominate the contributions originating from both the direct and Orbach pro- 
cesses. 

APPENDIX C: MAGNETOCRYSTALLINE ANISOTROPY 

In this appendix, we give a simple estimate of magnetocrystalline anisotropy in metals, 
which nevertheless contains most of the features of a complete bandstructure calculation 
of this quantity and already yields the correct order of magnitude. For that purpose, we 
consider a single, for simplicity parabolic, but spin-orbit split band (Fig. 4). Hereby we 
neglect the fact that parabolic bands usually represent s electrons which feel neither spin- 
orbit nor exchange interaction. In addition, we neglect the magnetic dipole-dipole coupling 
which favors in-plane magnetization in two dimensions and is zero in the bulk of cubic or 
hexagonal crystals such as Gd. It is in particular the spin-orbit induced magnetocrystalline 
anisotropy energy which may (but does not necessarily have to) favor a perpendicular easy 
axis in thin films and is therefore of interest for high-density magnetic recording (the time 
limit of which is related to tsl)- 

We calculate now the maximum energy gain from magnetocrystalline anisotropy in this 
model. This gain originates from the change of the band occupation up on spin-orbit induced 
lifting of the band degeneracy at the Fermi level. Electrons are transferred from one branch 
of the band to the other. Assuming a Brillouin sphere in three dimensions one therefore 
obtains an anisotropy energy of 



E 



'anisotropy 



= A 



s.o. 



X 




= A 



■s.o. 



X 



3AA; 



(3D) 



(28) 



where 



A A; 1 k 
kF 2 kp 



(29) 



12 



is the number of states contributing to the change of the electronic occupation. For typical 
values of the spin-orbit coupling constant A s o . = 70 meV and the Fermi energy E F = 10 eV 
we find the result 



E 



'anisotropy 



735/ieV (3D). 



(30) 



Thus, this crude model yields already important insights in some of the microscopic features 



model gives the correct order of magnitude for E an i SOtropy in films or non-cubic bulk crystals. 
The actual value for Gd might be somewhat smaller in Gd but this would even improve 
the agreement of t$l with experiment, (ii) The model shows that the magnetocystalline 
anisotropy is smaller than the spin-orbit coupling constant by two orders of magnitude since 
only a relatively small portion of states close to the Fermi level may gain energy from the spin- 
orbit induced lifting of degeneracies. For all other states the upward and downward shifts of 
the lifted degeneracies cancel, (iii) The model immediately yields that the anisotropy energy 
resulting from the lifting of degeneracies is proportional to the square of spin-orbit coupling, 
since besides the explicit linear dependence on A s o . also the portion of contributing states 
is linear in A s , .. There is no azimuthal dependence on spin-orbit coupling in this model 
in remarkable agreement with the line degeneracies found in Fe monolayers^. Thus, the 
energy gain resulting from the lifting of degeneracies close to the Fermi energy is of the 
same order of magnitude as the magnetocrystalline anisotropy originating from level shifts 
far below Ep, which can be obtained already in nondegenerate second-order perturbation 
theory with respect to spin-orbit coupling, (iv) The model explains why perpendicular 
anisotropy may be favored in thin films: Due to the reduced coordination number in these 
films, narrow bands of large density of states close to the Fermi level (in ferromagnets) may 
occur which can gain a suffciently large amount of magnetocrystalline anisotropy energy. 

(v) The anisotropy at interfaces and the nonlinear magneto-optical Kerr-Effect are closely 
related via spin-orbit coupling although the latter results from nonlinear optical excitations. 

(vi) Interface hybridization of a ferromagnet with a strong spin-orbit scatterer may yield 



of magnetocrystalline anisotropy which are confirmed 
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large anisotropies due to the reoccupation of many states close to their common Fermi level, 
(vii) Spin-orbit coupling does not split exchange-split bands again. Thus, the diagonal part 

L Z S Z (31) 

just yields contrary level shifts of spin- up and spin-down bands whereas the diagonal con- 
tribution 

L X S X + LySy = —(L+S— + L—S+) (32) 

yields spin-flips. It is these spin-flips that contribute to the spin-lattice relaxation time tsl 
which in this view describes the time required for adapting the (temperature dependent) 
magneto-crystalline anisotropy to the lattice temperature. Thus, our model yields the cor- 
rect order of magnitude and a change of the direction of the magnetic moments. 
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APPENDIX D: PHONON-MAGNON COUPLING 

If we want to conceive a similar theory for spin-lattice relaxation in ferromagnetic transition 
metals such as Fe we have to notice the following important difference: In the rare earth 
metal Gd it is sufficient to consider the localized /-electron spins carrying the magnetic 
moment of 7 hb, which gives rise to seven possible orientations {mi quantum numbers). 
Thus, Raman processes may easily take place. In the transition metal Fe, however, the d 
electrons are itinerant (delocalized), have to be described within the band picture, and have 
only two orientations (spin- up and spin-down). Therefore collective and quantized magnetic 
excitations (magnons) have to be allowed for the spin-lattice relaxation (Fig. 5). The 
transition Hamiltonian then describes the phonon-magnon coupling using boson creation 
and annihiliation operators aj^ for magnons and b-jp for phonons 

H ph _ mag (k) = Aj&kOk + h - c -) (33) 

withl 



^^f^^Ww (34) 

Here, D is some coupling strength, F some crystal field parameter, S the spin, M the effec- 
tive mass, and v s the averaged speed of sound. 

The interaction Hamiltonian H p h- m ag{k) can then be inserted in the full phonon-magnon 
Hamiltoniani 

H = Et<4«k + wE&t&k + H ph _ mag (k)] (35) 

k 

which is easily solved by applying the unitary transformation 

4° = AjpcosBk + B^sinGk 

&£) = B^cosGk - A£?srae k (36) 

to yield 
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uja = cos2 @k + c^k s i n2 ©k — 2D k cos 6k sin 6k 
u>a = cJ k a cos 2 6k + u; k sin 2 6k + 2Dk cos 6k sin 6k- 

6 k is given by 

tan 26k 



, ,P , ,m ' 



This then completes the formal solution of the phonon-magnon problem 
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FIGURES 

FIG. 1. (a) Direct process, (b) Orbach process, and (c) Raman process. 



FIG. 2. Microscopic model for magnetocrystalline anisotropy energy. 



FIG. 3. Spin-lattice relaxation in transition metals. 
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